Magnetic Helicity Conservation and Inverse Energy Cascade 
in Electron Magnetohydrodynamic Wave Packets 
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Electron magnetohydrodynamics (EMHD) provides a fluid-like description of small-scale magne- 
tized plasmas. An EMHD wave (also known as whistler wave) propagates along magnetic field lines. 
The direction of propagation can be either parallel or anti-parallel to the magnetic field lines. We 
numerically study propagation of 3-dimensional (3D) EMHD wave packets moving in one direction. 
We obtain two major results: 1. Unlike its magnetohydrodynamic (MHD) counterpart, an EMHD 
wave packet is dispersive. Because of this, EMHD wave packets traveling in one direction create 
opposite traveling wave packets via self-interaction and cascade energy to smaller scales. 2. EMHD 
wave packets traveling in one direction clearly exhibit inverse energy cascade. We find that the 
latter is due to conservation of magnetic helicity. We compare inverse energy cascade in 3D EMHD 
turbulence and 2-dimensional (2D) hydrodynamic turbulence. 



1. SIMULATIONS AND RESULTS 

Introduction — Astrophysical plasmas are observed 
in a wide range of length-scales. On large scales, we can 
treat such plasmas as conducting fluids and therefore 
we can use magnetohydrodynamics (MHD). Although 
MHD is a simple and powerful tool for large scales, it 
is not suitable for describing small-scale physics, espe- 
cially physics near and below the proton gyro-scale. Since 
many astrophysical processes critically depend on small- 
scale physics, proper description of small-scale physics 
is needed. There are several numerical models that can 
handle small-scale physics. Perhaps a full kinetic treat- 
ment would be the best for proper description of small- 
scale physics. However, a full kinetic description of plas- 
mas is still a challenge for modern computers. 

Electron magnetohydrodynamics (EMHD) is a fluid- 
like model of small-scale plasmas [lj and can be viewed 
as Hall MHD in the limit of kpi ^> 1, where pi is the 
ion gyroradius and k the wavenumber. On scales below 
the ion inertial length di = c/uj pi , where c is the speed of 
light and oj p i is the ion plasma frequency, we can assume 
that the ions create only smooth motionless background 
and fast electron flows carry all the current, so that 
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V x B, 



(1) 



where v e is the electron velocity, J is the electric current 
density, B is the magnetic field, n e is the electron number 
density, and e is the absolute value of the electric charge. 
Inserting this into the usual magnetic induction equation 
(dB/dt = V x (v e x B) + ??V 2 B), we obtain the EMHD 
equation 



dB 

~dt 
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V x [(V x B) x B] 



■7jV 2 B. 



(2) 



Note that, in this paper, we only consider the zero (nor- 
malized) electron inertial length case: d e — c/(uj pe L) — > 
0, where Lu pe is the electron plasma frequency and L is 
the typical size of the system. 



As is the case with large scales, magnetized turbu- 
lence on small scales also affects many physical processes 
and hence is of great interest for studies of magnetic 
reconnection[3, 3 space plasmas and the solar wind0- 
(l3j |. neutron stars fl4HL6j|. advection dominated accretion 
flows etc. Due to its simplicity, EMHD formalism has 



been used for studies of small-scale turbulence [18|- 25 1 . 
Earlier studies revealed that E b (k) oc fc~ 7 / 3 E IH l26j. 
and anisotropic turbulence structures (21, 23. 27|. 

In the presence of a strong mean field B , an MHD 
perturbation moves along magnetic field at the Alfven 
speed (oc Bq) and MHD wave packets moving in one di- 
rection do not interact each other and do not create tur- 
bulence. Therefore, collisions of opposite-traveling wave 
packets are essential for generation of MHD turbulence. 
In contrast, an EMHD perturbation moves along mag- 
netic field at a speed proportional to kB^ 1 which implies 
that a perturbation with larger k is faster than that with 
smaller k. As a result, whistler waves are dispersive and 
whistler wave packets moving in one direction can self- 
interact and produce small-scale structures (see for 
2D EMHD), which means that collisions of whistler wave 
packets are not essential for generation of EMHD turbu- 
lence. Therefore understanding the dynamics of EMHD 
wave packets is important for study of EMHD turbulence. 
Traveling EMHD wave packets can commonly occur in 
nature. Any local disturbances can create wave packets 
traveling along magnetic field lines. Therefore, propa- 
gation of whistler wave packets moving in one direction 
deserves a scrutiny. In this paper, we study propagation 
of 3D EMHD wave packets in detail. 

Numerical Method — We have calculated the time 
evolution of 3D incompressible EMHD wave packets mov- 
ing in one direction. We have adopted a pseudospectral 
code to solve the normalized incompressible EMHD equa- 
tion in a periodic box of size 2tt: 

-=-Vx[(VxB)xB]+,'V 2 B, (3) 
at 

where the magnetic field, time, and length are normalized 
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by a mean field Bq, the whistler time t w = L 2 (ui pe / c) 2 /O e 
(fi e = electron gyrofrequency) , and a characteristic length 
scale L (see, e.g., (Hj). The resistivity rj' in equation ((3]) 
is dimensionless. The dispersion relation of a whistler 
wave in this normalized unit is u) — kk\\Bo, where k\\ 
is the wave number parallel to the mean magnetic field. 
The magnetic field consists of the uniform background 
field and a fluctuating field: B = Bo + b. The strength 
of the uniform background field, Bq, is set to 1. We use 
either 256 3 or 512 3 collocation points. At t = 0, all waves 
are moving in the same direction and their wave numbers 
are restricted to the range 8 < fc < 15 in wavevector 
(k) space. The direction of propagation corresponds to 
the positive direction of the mean magnetic field in our 
simulations. (Hereinafter, we use positive and negative 
to denote the direction of wave propagation with respect 
to the magnetic field.) The amplitudes of the random 
magnetic field at t — is ~ 1. Hyperdiffusivity is used for 
the diffusion term. The power of hyperdiffusivity is set 
to 3, so that the dissipation term in the above equation 
is replaced with 77 3 (V 2 ) 3 B, where 773 is approximately 
3 x 10" 10 for 256 3 and 1 x 1CT 11 for 512 3 . 

Results — Figure Q] shows time evolution of magnetic 
energy density. The solid curves denote the total mag- 
netic energy density (i.e. the total energy of the wave 
packets moving in both the positive and the negative di- 
rections) and the dotted curves the energy of wave pack- 
ets moving in the negative direction. We can see that 
wave packets initially moving only in the positive direc- 
tion can create waves moving in the negative direction. 
Although the energy of wave packets moving in the neg- 
ative direction is small (at most a few % compared with 
the energy of the wave packets moving in the positive 
direction), we clearly observe that magnetic energy de- 
cays. Results from 512 3 (thick curves) and 256 3 (thin 
cures) show a reasonable agreement. 

Figure [2] shows magnetic energy spectrum as a function 
of time. At t = 0, Fourier modes between k = 8 and k = 
15 are excited. The long-dashed curves in Figure [2] show 
the initial spectrum. As time goes on, the initial energy 
cascades down to smaller scales and, as a result, a power- 
law-like spectrum forms for k > 15. At the same time 
the peak of the energy spectrum moves to larger scales, 
so that the wavenumber at which the spectrum peaks, 
k p , gets smaller. We clearly observe inverse cascade of 
magnetic energy. 

The magnetic helicity is a conserved quantity in 
EMHD 28;]. Figure [3] shows that the magnetic helicity 
is extremely well conserved. The solid curves denote the 
net magnetic helicity (i.e. the helicity of wave packets 
moving in the positive direction minus that in the nega- 
tive direction) and the dotted curves the helicity of wave 
packets moving in the negative direction. The helicity 
of the wave packets moving in the negative direction is 
much smaller than that of the waves moving the other 
way. As in Figure [TJ thick curves are for 512 3 and thin 



curves for 256 3 . 

2. THEORETICAL INTERPRETATIONS 

Helicity decomposition — An EMHD wave travel- 
ing along the magnetic field lines is circularly polarized. 
In Fourier space, the bases that describe circular polar- 
ization are 

e+ = (§1 + is 2 )/\ / 2 and e_ = (§1 - is 2 )/V2, (4) 

where §1 and S2 are orthogonal unit vectors perpendicu- 
lar to k and we assume siXS2 = k = k//c. Any EMHD 
Fourier mode can be decomposed into a '+' wave and a 
'-' wave: 

b k = M+ + b-i-, (5) 

where bk = bi&i + 62S2, b+ = (b\ — ibi)j\J 1 l 1 and £>_ = 
(61 + fb 2 )/V2. 

The '+' wave moves in the positive direction and the '-' 
wave in the negative (i.e. apposite) direction with respect 
to magnetic field. Let us consider a '+' wave. In Fourier 
space, we can show that 

d(b+e+)/dt = -ikk\\B (b+i+) (6) 

(see, e.g. [11]). Therefore, we have b + = b exp(— iut), 
where bo is a constant. Fourier transform of 6+ to real 
space gives a plane wave whose velocity relative to the 
magnetic field is positive. 

The '+' wave has a positive magnetic helicity and the 
'-' wave a negative magnetic helicity. Let us consider a 
'+' wave. In Fourier space, we can show that the Fourier 
component of the vector potential is given by 

a k = a + e + = (b + /k)e + , (7) 

where we use the Coulomb gauge. Therefore, the mag- 
netic helicity of the '+' wave is positive. 

Inverse cascade — The magnetic helicity is a con- 
served quantity in EMHD (see, e.g. [HI]). For an EMHD 
wave packet composed of only '+' waves, the magnetic 
helicity spectrum Eh{k) is simply 

E h {k) = E h {k)/k. (8) 

Therefore, the helicity dissipation rate becomes negligi- 
ble when k ( i, where kd is the dissipation wavenumber, 
is large with respect to k p . Note, however, that energy 
dissipation rate can be non-negligible. Therefore, as en- 
ergy dissipates, the peak wavenumber k p should become 
smaller. 

Since J Eh(k)dk is nearly constant and J Eh(k)dk ~ 
[Eb(k p ) I kp]k p ~ Eb(k p ), we have Eb(k p ) ^constant (see 
Figure 2; see also left panel of Figure HJ. Since b 2 ~ 
Eb(k p )kp oc kp, we expect that 

k p {t) cx b 2 {t), (9) 
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FIG. 1. Time evolution of the fluctuating magnetic energy density b 2 . Initially all wave packets move in the same direction. 
Due to self-interaction, b 2 decreases as time goes on (solid curves). Self-interaction generates waves traveling in the opposite 
direction (dotted curves). The energy of the opposite-traveling waves is no larger than a few % of the energy of the waves 
moving in the original direction. The thin curves are for 256 3 and the thick ones for 512 3 . 



FIG. 2. Magnetic spectra showing inverse energy cascade. The peak of magnetic energy spectrum moves to larger scales. Note 
that E(k p ), where k p is the wave number at which the magnetic energy spectrum peaks, is almost constant. Magnetic helicity 
conservation plays a key role in the inverse cascade. 
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FIG. 3. Helicity conservation. The net magnetic helicity 
(solid curves) is well- conserved. The dotted curves denote 
the absolute values of magnetic helicity in opposite-traveling 
wave packets, which are generated by self-interaction of the 
initial wave packets. The thin curves are for 256 3 and the 
thick ones for 512 3 . 



which is confirmed by right panel of Figure 2] We ob- 



tained k p in such a way that X)fc-o.5<|k'|<fc+o.5 l D (k')| 2 
(k = 1.0, 1.2, 1.4, 1.6, ...) is maximum at k = k p . 

Comparison with 2D hydrodynamic turbulence 

— We find a very good correspondence between 3D 
EMHD turbulence and 2D incompressible hydrodynamic 
turbulence. The 2D incompressible hydrodynamic equa- 
tion has two ideal invariants: the energy and the en- 
stropy (oc J |V x v| 2 <i 2 x). The spectrum of enstropy is 



E h (k) 




FIG. 4. Quantities that remain constant during inverse en- 
ergy cascade. We expect that Eb(k p ) and b 2 (t)/k p {i) are al- 
most constant. Left panel: Schematic diagram showing con- 
stancy of Ef>(kp). See Figure [5] for actual simulation results. 
Right panel: Simulation result that confirms constancy of 
b 2 (t)/k p (t). 



E en {k) — k 2 E v (k), where E v (k) is energy spectrum of 
velocity. In driven 2D hydrodynamic turbulence, both 
energy and enstropy involve turbulence cascade: en- 
stropy exhibits forward cascade and energy inverse cas- 
cade. When energy is the cascading quantity, which is 
the case for scales larger than the energy injection scale, 
we have 



kE v (k) 



I/Vt l/(ky/kEJk) 



=const. — > E v (k) oc k 



-5/3 



(10) 



where we used vi ~ y/kE v (k). On the other hand, when 
enstropy is the cascading quantity, which is the case for 
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scales smaller than the energy injection scale, we have 
kk 2 E v (k) 



=const. — !• Ey(k) oc k 



(11) 



l/(k^/kE v {k) 

where we used E en (k) = k 2 E v (k). In this case, interac- 
tions are non-local in wavevector space. 

In 3D EMHD turbulence, magnetic energy and mag- 
netic helicity are conserved quantities in the absence of 
dissipation. Therefore, either magnetic energy or mag- 
netic helicity can involve energy cascade. Using argu- 
ments similar to the 2D hydrodynamic case, we can ob- 
tain magnetic energy spectrum. When magnetic energy 
is the cascading quantity, we have 

kE h (k) 



l/vi l/(k^/kk 2 E b (k) 



=const. 



E b (k) oc £r 7 / 3 , 



(12) 

where we used v oc J oc V x B, which means E v (k) oc 
k 2 E b (k). On the other hand, when magnetic helicity is 
the cascading quantity, we have 

kE b (k)/k 



=const. 



E b (k) oc fc" 5/3 , (13) 



l/(k^kk 2 E b (k) 

where we used Eq. ([5]). 

When we inject magnetic helicity (and, hence, mag- 
netic energy) on a scale, we expect to see both inverse 
and forward energy cascade. It is well-known that when 
energy is injected on a scale, the energy cascades down to 
smaller scales and the small-scale magnetic energy spec- 
trum is proportional to A; -7 / 3 . Therefore, it is evident 
that magnetic energy exhibits forward cascade. However, 
it is not clear what is the cascading entity for inverse cas- 
cade. We will address this issue elsewhere (Kim et al., in 
preparation) . 

In summary, we have shown that EMHD wave pack- 
ets moving in one direction can create opposite-traveling 
wave packets through self-interaction and that, because 
of magnetic helicity conservation, EMHD wave packets 
moving in one direction show inverse energy cascade. In- 
verse cascade of energy can affect transport phenomena, 
such as heat transport, and potentially be a source of 
plasma instabilities. It is also potentially important for 
magnetic reconnection and evolution of magnetic field in 
neutron stars. In general, traveling EMHD wave packets 
can be a source of magnetic helicity, which could affect 
evolution of a large-scale magnetic field. 

This research was supported by National R&D Pro- 
gram through the National Research Foundation of Ko- 
rea(NRF) funded by the Ministry of Education, Science 
and Technology (No. 2010-0029102). JC thanks H. Kim 
for useful discussions. 



jcho@cnu.ac.kr; Dept. of Astronomy & Space Science, 
Chungnam National Univ., Daejeon, Korea 



[8 

[9 

[10 
[11 

[12 
[13 
[14 

[is; 

[16 

[17; 

[18 

[19 

[20 

[21 

[22 

[23 
[24 

[25 

[26 
[27 
[28 



A. S. Kingsep, K. V. Chukbar, and V. V. Yan'kov, Re- 
views of Plasma Physics, Vol. 16 (Consultants Bureau, 
New York, 1990). 

J. F. Drake, R. G. Kleva, and M. E. Mandt, 

physical Review Letters 73, 1251 (1994)| 

D. Biskamp, |Astrophys Space Sci. 242, 165 (19 96) 

P. Dmitruk~ and W. H Matthaeus, 



and H. 



Li. 



Physics of Plasmas 13, 042307 (2006) 
S. P. Gary, S. Saito, 
|Geophys. Re. Lett. 35, 2104 (2008)| 
G. G. Howes, W. Dorland, S. C. Cowley, G. W. 
Hammett, E. Quataert, A. A. Schekochihin, and T. Tat- 
suno, physical Review Letters 100, 065004 (2008) | 
larXiv:0711.4355l 

G. G. Howes, S. C. Cowley, W. Dorland, G. W. Ham- 
mett, E. Quataert, A. A. Schekochihin, and T. Tatsuno, 
physical Review Letters 101, 149502 (2008) | 
R. J. Leamon, C. W. Smith, N. F. Ness, 

W. H. Matthaeus, and H. K. Wong, 

|J. Geophys. Res. 103, 4775 (1998)| 

W. H. Matthaeus, S. Servidio, and P. Dmitruk, 



Physical Review Letters 101, 149501 (2008) 



arXiv:0812.2493. 

S. Saito, S. P. Gary, H. Li, and Y. Narita, 

physics of Plasmas 15, 102305 (2008)| 

A. A. Schekochihin, S. C. Cowley, W. Dorland, G. W. 

Hammett, G. G. Howes, E. Quataert, and T. Tatsuno, 

|Astrophys. J. Supp 182, 310 (2009)| larXiv:0704.0044l 

C. W. Smith, K. Hamilton, B. J. Vasquez, and R. J. 

Leamon, |Astrophys. J. Lett. 645, L85 (2006)] 

O. Stawicki, ST P Gary, and H. Li, 

| J. Geophys. Res. 106, 8273 (2001) | 

P. Goldreich and A. Reisenegger, 

|Astrophys. J. 395, 250 (1992) | 

A. Cumming, P. Arras, and E. Zweibel, 

|Astrophys. J. 609, 999 (2004)| |arXiv:astro-ph/0402392| 
C. J. 



Wareing and 



R. 



Astron. Astrophys. 508, L39 (2009) 



arXiv:0912.2656 [astro-ph.HE] 
Quataert 



E 



and 



A. 



Hollerbach, 



Gruzinov, 



|Astrophys. J. 520, 248 (1999) | |arXiv:as tro-ph/98 63112[ 
D! Biskamp, W. Schwarz, and J. F. Drake, 
physical Review Letters 76, 1264 (1996)| 
D. Biskamp, E. Schwarz, A. Zeiler, A. Celani, and J. F. 
Drake, physics of Plasmas 6, 751 (1999)| 

C. S. Ng, A. Bhattacharjee, K. Germaschewski, and 
S. Galtier, pPhy"sics of Plasmas 10, 1954 (2003) | 

J. Oho and ~K. Lazarian, 

|Astrophys. J. Lett. 615, L41 (2004)] 

D. Shaikh and P. Zank, 

physics of Plasmas 12, 122310 (20"05)| 

J. Cho and A. Lazarian, |Astrophys. J. 701, 236 (2009)| 



Wareing and 



R. 



Physics of Plasmas 16, 042307 (2009) , 
arXiv:0904.4825 [physics. plasm-ph] 
C. T. Wareing and R. 



Hollerbach, 



Hollerbach, 



Journal of Plasma Physics 76, 117 (2010) 
arXiv:0906.3401 [physics. plasm-ph] 

S. I. Vainshtein, Sov. Phys. JETP 37, 73 (1973), [Zh. 
Eksp. Teor. Fiz. 64, 139 (1973)]. 

S. Dastgeer, A. Das, P. Kaw, and P. H. Diamond, 
Physics of Plasmas 7, 571 (2000) | 

5! Galtier and A. Bhattacharjee, 

Physics of Plasmas 10, 3065 (2003)] 



